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Apology 

While this article is dedicated to both Mourad Ismail and Dennis Stanton, it does not directly 
reference any of their works. The main reason is that I only talk about the most trivial kind of 



OO 



recurrences: linear and constant coefficients. But try searching Ismail AND Recurrences or 
Stanton AND Recurrences in the database MathSciNet (or Google Scholar) and you would see 
that both Mourad and Dennis are great gurus in recurrences, so the subject matter of this paper 
is not entirely inappropriate as a tribute to them. The present work is also largely experimental, 
and Dennis Stanton is a great pioneer in computer experimentations! 
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jA ! PROLOGUE 

a . 

Before starting the paper itself, let me very briefly mention what I talked about at the wonderful 
conference 

^. ■ q-Series 2011: An International Conference on q-Series, Partitions and Special Functions, 

m 

"honouring Mourad Ismail and Dennis Stanton for their valuable contributions to Number 
theory and Special Functions throughout their careers", that took place on March 14-16, 2011 
at Georgia Southern University, and perfectly organized by Drew Sills. My plenary talk was 



o 



% 



entitled 

u Some Golden Oldies of Mourad Ismail and Dennis Stanton" 

and the abstract was short and sweet: 

"This is the time to recall some of the beautiful mathematics that I learned from Dennis and 
Mourad". 

I talked, inter alia, about Dennis Stanton's amazing article 

"A short proof of a generating function for Jacobi polynomials" , 
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that appeared in Proc. Amer. Math. Soc. 80 (1980), 398-400, where he used an idea that he 
attributed to his advisor, Dick Askey, but that Askey modestly claims goes back to Hermite, to 
give a very soft and elegant proof of a very hairy formula of Bailey 

Another thing that I talked about was Mourad Ismail's gorgeous article 

"A simple proof of Ramanujan's i^i sum" , 

that contained the proof- from-the-book of Ramanujan's lovely formula 
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(where, as usual (A) n := (1 - A) ■ ■ ■ (1 - q n ~ 1 A), if n > 0, (A) n := 1/((1 - g"M) • • • (1 - g n A)), if 
n < 0, and (A)^ := (1 - A)(l - qA)(l - q 2 A) ■■■) . 

Mourad's article appeared in the proceedings (of the AMS) three years earlier (PAMS 63 (1977), 
185-186.) This proof was already immortalized and canonized (and shrunk to half a page!) in 
Appendix C of George Andrews' classic monograph, CBMS #66, that appeared in 1986 and was 
based on ten beautiful lectures delivered at an NSF-CBMS conference that took place in 1985, at 
Arizona State University, and organized by Mourad Ismail and Ed Ihring. 

In my talk, I mentioned that while Mourad's proof certainly qualifies to be included in "God's 
book" , since, like Mourad and George, God is, by definition, an infinitarian, it does not qualify to 
be included in my book. It said (in George's rendition, my emphasis) 

"Regarding the left-hand side as an analytic function of b for |6| < 1 . . . " . 

Then one plugs-in b = q N (TV = 0, 1, 2, . . .), getting the trivial q-binomial theorem, and one sees 
that the left side minus the right side vanishes for "infinitely" many values of b and then uses the 
"fact" that an "analytic" function inside |6| < 1 that vanishes on a "convergent" "infinite" sequence 
"must" be identically zero. 

Of course, to finitists like myself, this proof is entirely non-rigorous, since it uses fictional things 
like so-called analytic functions, and uses heavy guns from a sophisticated (and flawed!) "infinite" 
theory. 

But don't despair! It is very easy to translate Mourad's flawed proof and make it entirely legit, 
and in the process make it even nicer. Replace the phrase 

"analytic function of b defined in \b\ < 1" 

by 

"bilateral formal power series in t whose coefficients are rational functions of b" , 
and note that the difference of the left and right sides is a bilateral formal power series (in t) 
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whose coefficients are rational functions of b (and of course also of a and q but that's irrelevant). 
A rational function of b, whose degree of the numerator is, say, m, is identically zero if it vanishes 
at m + 1 distinct values of 6, so the "infinitely" many points b = q N (N = 0, 1, 2 . . .) is more than 
enough. 

FEATURE PRESENTATION 

The Maple package Cf inite 

This article is accompanied by the Maple package Cf inite downloadable directly from 

http : //www .math. rutgers . edu/~zeilberg/tokhniot/Cf inite, 

or from the "front" 

http://uww.math. rutgers . edu/~zeilberg/mamarim/mamarimhtml/cf inite.html, 

where one can find links to fiftenn sample input and output files, some of which are mentioned 
throughout this paper. 

C-finite Sequences 

Recall that a C-finite sequence {a(n)}, n = 0, 1, ... is a sequence that satisfies a linear-recurrence 
equation with constant coefficients. It is known (but not as well-known as it should be!) and 
easy to see (e.g. [Z2],[KP]) that the set of C-finite sequences is an algebra. Even though a C- 
finite sequence is an "infinite" sequence, it is in fact, like everything else in mathematics 
(and elsewhere!) a finite object. An order-L C-finite sequence a{n) is completely specified by the 
coefficients c\ , Ci , . . . , Cl of the recurrence 

a(n) = c\a{n — 1) + C2<a(n — 2) + . . . + CLa(n — L) , 

and the initial conditions 

a(0) = d\ , ... , a(L — 1) = cIl 

So a C-finite sequence can be coded in terms of the 2L "bits" of information 

[[di,...,d L ],[ci,...,c L ]] ■ 

For example, the Fibonacci sequence is written: 

[[0,1], [1,1]]. 

Since this ansatz (see [Z2]) is fully decidable, it is possible to decide equality, and evaluate ab initio, 
wide classes of sums, and things are easier than the holonomic ansatz[Zl]. The wonderful new book 
by Manuel Kauers and Peter Paule[KP] also presents a convincing case. See [HX] [GW] [Kau] [KZ] 
for very interesting and efficient algorithms. 



Rational Generating Functions 

Equivalently, a C-finite sequence is a sequence {a(n)} whose ordinary generating function, 
X^^Lo a ( n ) z " '> * s a ra tional function of z. These come up a lot in combinatorics and elsewhere 
(e.g. formal languages). See the old testament[St2], chapter 4, and the new testament[KP] , chapter 

4. 

[To go from a C-finite representation to a rational function, use CtoR(C,z) ; in Cf inite. To go 
the other way, do RtoC(f ,z) ; .] 

Etymology 

I coined the term C-finite sequence in [Zl], as a hybrid analog of Richard Stanely's [Stl] names "In- 
finite function" and "P-recursive sequence" . If I had to do it over I would call them "C -recursive 
sequences" , but it is too late now since the term C-finite already made it into the wonderful 
undergraduate textbook [KP], and it is also in the title of the important paper [GW]. 

Zeilberger-style proofs: You (Often) CAN generalize from FINITELY Many Cases 

The conventional wisdom of mathematics (at least for the last 2500 years), preached to us by our 
teachers and that, in turn, we preach to our students, is that you can't generalize from finitely many 
cases. While this is certainly true sometimes, it is not always true. Many times you can generalize 
from finitely many cases, just like natural scientists. 

Michael Hirschhorn kindly called this style of proof "in the spirit of Zeilberger", see his beautiful 
proof[H] of an amazing identity of Ramanujan, that gives infinitely many "almost" counterexamples 
to Fermat's Last Theorem for n = 3, namely infinitely many triples {(a, b, c)} such that a 3 +b 3 +c 3 = 
±1. 

But not Everyone Knows About this Style of Proof 

Everybody knows that numerical identities like 2 + 7 = 3x3 are routinely provable, using standard 
algorithms. But even people like Neil Sloane (and Jeff Shallit) and James Sellers are not fully aware 
that identities amongst C-finite sequences are equally routine. If they were, James Sellers' article 
[Se] would not have been accepted for publication in Journal of Integer Sequences, or written in 
the first place, see the parody [Sr]. There are hundreds (possibly thousands) of articles like this in 
the literature, sometimes giving "elegant" proofs of such trivial results. While it is always nice to 
have elegant proofs, honesty requires that the authors state clearly, in the abstract, that the result 
that they are elegantly proving is routinely provable. See my opinion [Z4]. 

The First Reason for This Article: Educating 

Since even Neil Sloane, Jeff Shallit, James Sellers, Art Benjamin, and many other people are not 
(fully) aware of the triviality of the C-finite ansatz (or more politely, there being an algorithmic 
proof theory for it), and in spite of the articles and book cited above, I thought that it is a good 



idea to make it better known. 

The Second Reason for This Article: Implementation 

While Curtis Greene and Herb Wilf [GW] and Manuel Kauers[Kau] (see also [KZ]) already have 
Mathematica implementations of many operations on C-finite sequences, and possibly also Maple 
ones, I thought that it is a good idea to design a C-finite calculator, that also enables one to 
discover new identities. The novelty, in that part, is the approach, pure guessing! (that is justified 
a posteriori). 

The Third Reason for This Article: Factorization 

The truly novel part (I believe) is in addressing the problem of factorization. See below. 

The C-finite Calculator 

In order to decipher a C-finite sequence where the first few terms are given, all you need is use linear 
algebra to "guess" (using the ansatz) the c's (you already know the d's). See Procedure GuessRec 
in Cf inite. If you have two C-finite sequences C\ and Ci of order L\ and L 2 , you don't need any 
fancy footwork to figure out the sequence C\ + C2 (of order L\ + L 2 ) and the sequence C1C2 (of 
order L\L 2 ). All you need is to crank out L\ + L 2 + 4 (the +4 is for safety reasons) and L\L 2 + 4 
terms, respectively, and let the computer guess the C-finite description, completely by guessing, 
using undetermined coefficients that is implemented by Procedure GuessRec in Cf inite. 

So it is (very!) easy to multiply C-finite sequences, in other words, from the C-finite descriptions 
of C\ and C 2 get the C-finite description of C1C2 (by C1C2 we mean the sequence whose n-th term 
is Ci(n)C2(n), in terms of their generating functions it is called the Hadamard product). 

Procedures Khibur and Kef el of Cf inite 

In order to add two C-finite sequences CI and C2, simply type 

Khibur (CI ,C2); 

In order to multiply type: 

Kef el (CI, C2); 

and for a verbose version, presenting fully detailed proofs in the spirit of Zeilberger, type 
KefelV(Cl,C2); . 

An Example of Using Procedure Kef el 

Let {U n (x)}^ =0 be the sequence of Chebyshev polynomials of the second kind, i.e. the sequence of 



polynomials in x defined by 

*-" y ' l-2xt + t 2 

n=0 

Typing 

CtoR(Kef el(Ux(a) ,Ux(b)) ,t) ; would give you, in a few seconds, the following result[Sh] 

Lou Shapiro's Product-Of-Two-Chebyshev-Polynomials Identity 

l-t 2 



Y,U n (a)U n (b)f 



l-4a6t-(-4a 2 + 2-46 2 )t 2 -iabt 3 + t A 

n=0 

Typing 

CtoR(Kefel(Kefel(Ux(a) ,Ux(b)) ,Ux(c)) ,t) ; 

yields, in a few more seconds, the following much deeper result 

Shalosh B. Ekhad's Product-Of-Three-Chebyshev-Polynomials Identity 

f; U n (a)U n (b)U n (c)t n = ^ , 

where the polynomials N(t) are D(t) are as follows. 

N(t) = I + (-4a 2 - 4b 2 - 4c 2 + 3) t 2 + IQ abet 3 + (-4a 2 - 46 2 - 4c 2 + 3) t 4 + t 6 , 
D(t) = t 8 -8 abet 7 + (16 a 2 b 2 + 16 a 2 c 2 -8a 2 + 16b 2 c 2 -8b 2 -8c 2 + A) t 6 + (-32 a 3 be + 40 abc - 32 ab 3 c - 32 abc 3 ) 
(16 a 4 + 64 a 2 b 2 c 2 - 16 a 2 + 16 o 4 - 16 b 2 + 6 + 16 c 4 - 16 c 2 ) t 4 + (-32 a 3 bc + 40 a6c - 32 ab 3 c - 32 abc 3 ) t 3 
+ (16 a 2 fe 2 + 16 aV - 8 a 2 + 16 6 2 c 2 - 8 6 2 - 8 c 2 + 4) t 2 - 8 abet + 1 . 

The (Computationally) Hard Problem of Factoring C-finite Sequences 

Alas, going backwards (just like in integer factorization, that makes our ATM cards hopefully 
secure) is not so easy! If you are given a C-finite sequence of order L, say, and L is composite, 
L = L1L2, (with Li,L,2 > 1) you would like to know whether there exist C-finite sequence C\ 
and C2 such that C = C1C2, and if they do, find them. One way, that works for small L, is to do 
symbolic multiplication of generic C-finite sequences, and then try to solve, by matching coefficients, 
the resulting non-linear system of algebraic equations, using the Buchberger algorithm. [This is 
implemented in procedure Factor ize of the Maple package Cf inite .] But for larger orders this 
is hopeless! 
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Procedure Factorizell does the same by brute force, but only handles integer sequences. While 
it can't go very far, it discovered, ab initio, in less than half a second, the three factorizations in 
[Se]. See 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/oCf inite4 

For a more verbose version see 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/oCf inite5 , 

and for a completely spelled-out proofs, in the spirit of Zeilberger, see 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/oCf inite5a 

Why is this problem interesting? 

Two great landmarks of Statistical Physics are the Onsager[0] solution of the two-dimensional 
Ising model and the Kasteleyn[Kas]- Temperley-Fisher[TF] solutions of the dimer problem. They 
use lots of human (ad-hoc) ingenuity to first get an explicit answer for a finite strip of arbitrary 
(symbolic) width. They then take the so-called thermodynamic limit. It turns out that in either 
case, the m-wide strip sequence is a C-finite sequence of order 2 m . Surprisingly, in both cases they 
happen to be products of m C-finite sequences of order 2 (different, but closely related). 

Since nowadays computers can automatically, completely rigorously, figure out the C-finite 
description for each strip of width m, for specific, numeric m, ([EZ],[Z3]) (in practice easily for 
m < 10), knowing how to "factorize" them explicitly, would lead one to conjecture both solutions, 
with fully rigorous proofs for m < 10, and with larger computers, beyond. Since physicists are not 
as hung-up as mathematicians about rigorous proofs, that would have been a great breakthrough, 
even without the human proofs for general m. Besides, the explicit "conjecture" discovered by the 
computer might suggest and inspire (to obtuse mathematicians) a formal proof. 

The "Cheating Algorithm" 

Since it is so hard to factorize explicitly, it is still nice to know, as fast as possible, whether or not 
the inputted C-finite sequence is factorizable. If it is not, it would be stupid to waste efforts in 
trying to factorize it. If it is, then it is worthwhile applying for time on a bigger computer. 

Recall (Binet) that "generically" every C-finite sequence [[d\, . . . , di], [c\, . . . , cl]] of order L can 
be written as a linear combination 



C(n) = Y,C, 



a ; 



1=1 



where the Cj's depend on the initial conditions, and the ckj's are the roots of the characteristic 
equation 



L 
L 



Z 

i=l 



J>* L-i = o 
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So if C := C1C2 and the roots of C\ and C 2 are ai, . . . , a/ /1 , and /?i, . . . , /3^ 2 respectively, then the 
roots of C, let's call them 71, . . . ,^LiL 2 i consist of the Cartesian product 

{ai&|l<i<Li , l<j<£ 2 } • 
with L := L1L2 elements. 
If this is indeed the case, then the set of I? ratios 

{-\l<i,j<L} , 

li 

would have a certain profile 0/ repetitions that the computer can easily figure out for arbitrary 
symbols ct\ , . . . , a^ x and f5\ , . . . , /3l 2 . 

Procedure Prodlndictor of the Maple package Cf inite 

Procedure ProdIndicator(m,n) yields the profile of repetitions indicative of the characteristic 
roots of a C-finite sequence that happens to be the product of a C-finite sequence of order m and 
a C-finite sequence of order n. 

For example, ProdIndicator(2,2) yields: 

[1,1,1,1,2,2,2,2,4] . 
To understand what is going on, let's work it out by hand. 

7i = «i/5i , 72 = ai#2 » 73 = a 2 Pi , 74 = a 2 (3 2 . 



So 



7i "i/5i 1 7i oil Pi Pi 7i oi/5i «i 71 «i/3i 



7i «i/?i 72 ai/?2 /?2 73 Oi 2 P\ a 2 74 «2/?2 

72 = Qi/^2 = & 72 = Q1A2 = 1 72 = Qi/^2 72 = Qiffa = «i 
71 ai/?i /?i 72 ai/3 2 73 a 2 /3i 74 «2/?2 «2 

73 _ a 2 Pi _a 2 73 _ a 2 /3i 73 _ a 2 Pi _ 73 _ «2/3i _ Pi_ 
71 «i/?i «i 7 2 ai(3 2 73 a 2 /?i 74 a 2 f3 2 Pi 

74 02/?2 74 ct 2 /3 2 a 2 74 a 2 /3 2 /3 2 74 a 2 j3 2 



7i "i/?i 72 "i/?2 «i ' 73 a 2 /3i /?i ' 74 a 2 /3 2 
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We see that the multi-set of all 16 ratios has: four occurrences of 1, two occurrences each of 4 1 ,4 2 -, 

' pi pi 

37,37; and four singletons, namely ° g , ° 2 ^ and their reciprocals. 

Now for the proposed C-finite sequence of order L, find (in floating point!, but with Digits : =100 ; ) 
approximations to the roots of its characteristic equation, then form these L 2 ratios, and group 
them into classes with the "same" value (up to the agreed-on approximation) . If you get the same 
pattern of repetition, then you have proved (empirically) that the given C-finite sequence C, of order 



L = L1L2, is indeed the product of C-finite sequences of orders L\ and L 2 . Procedure IsProd in 
Cf inite implements this algorithm. See the source code for more details. 

If C has order L = L 1 L 2 • • • L r , and you want to find out whether C is a product of r C-finite 
sequences of orders Li, . . . , L r you do the analogous thing. Procedure IsProdG in Cf inite imple- 
ments this more general scenario. 

Output 

Using the output from [EZ] we confirmed that the straight enumeration dimer problems for strips 
of width < 10 are indeed products of C-finite sequences of order 2. See 
http : //www . math . rutgers . edu/~zeilberg/tokhniot/oCf inite2 

Using the output from [Z3] we confirmed that the weighted enumeration dimer problem for strips 
of width < 10 indeed are products of C-finite sequences of order 2 for many random numerical 
assignments of the weights. 

As for the actual factorization, we were, on our modest computer, only able to find them (from 
scratch, without peeking at the answer) for m < 7, but a more clever implementation, and a larger 
computer, no doubt would be able to conjecture ab initio (without any human ad-hocery!) the 
exact solution of the dimer problem derived and proved in [Kas] and [TF]. Ditto for Onsager's [O] 
(human) tour- de- force. 

Sample output for some of the other procedures (e.g BT for the Binomial Transform and GuessNLR 
for finding non-linear (polynomial) recurrences of lower-order than the (linear) order of a given 
C-finite sequence) can be obtained from the front of the present article: 

http : //www . math . rutgers . edu/~zeilberg/mamarim/mamarimhtml/cf inite . html . 

Encore: 1142 beautiful and deep Greene- Wilf-style Fibonacci identities in less than 
4400 seconds 

See the computer-generated webbook 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/oCf initel3 

For simpler identities see also 

http: //www. math. rut gers.edu/~zeilberg/tokhniot/oCf inite 11 and 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/oCf initel2 . 

These fascinating new identities can keep bijective combinatorialists busy for the next one hundred 
years. Each of these identities cries out for an insightful, elegant, combinatorial proof! 
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